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A THEORETICAL AND EXPERIMENTAL STUDY OP HYPERSONIC FLOW 


OVER FLARED BODIES AT INCIDENCE 

By John V. Rakich and Gene P. Menees 
Ames Research Center 


SUMMARY 


The inviscid flow over flared bodies at angle of attack is investigated. 
Experimental pressure distributions are presented for bodies with a hemispheri- 
cal nose and with a conical nose. Data are shown for Mach numbers of 5 *2, 7*^ 
and 10.5 and for angles of attack of 0°, 5°, and 10°. The models are provided 
with a means for removing the boundary layer upstream of the flare to eliminate 
shock-induced separation. 

A perturbation method is developed for the numerical calculation of the 
flare shock conditions, and the resulting equations are incorporated into a 
computer program based on the linearized characteristics method. Comparisons 
are made with experiment and with other more approximate theories which have 
been developed for flared bodies. 


INTRODUCTION 


Aircraft and spacecraft designers are often faced with the problem of 
predicting the aerodynamic forces on bodies that have compression corners on 
the surface. A particular example is the class of flare -stabilized entry 
bodies made up of a cylinder followed by a conical frustum. If the flow 
approaching the flare is locally supersonic, a shock wave is formed near the 
corner between the cylinder and flare. This shock interacts with the boundary 
layer on the cylinder ahead of the corner, and the boundary layer perturbs the 
flow approaching the flare. Thus, calculation of the inviscid and boundary 
layer flow is coupled, making the analysis of both more difficult. The embed- 
ded shock can cause the boundary layer to separate, and can therefore have a 
large effect on the inviscid flow. The aerodynamic stability at small angles 
of attack is especially sensitive to this interaction. Even the development 
of uncoupled Inviscid theories has been hindered by the fact that experimental 
results have usually had various degrees of boundary-layer shock -wave 
interaction. 

Numerous attempts have been made to predict the inviscid flow over flared 
bodies at incidence (see, e.g., refs. 1 to 5). These existing theories have 
tended to be empirical and have consisted of patching together various simple 
flows. While adequate for most preliminary design purposes, these theories 
have some deficiencies which appear in comparisons with experiment (see ref. 2), 
and it is not possible to determine whether the differences are due to approx- 
imations in the theory or to viscous effects in the experiments. The present 
paper therefore has a dual purpose. First, an attempt is made to apply exact 
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numerical methods to determine the perturbations to the flow over flared 
bodies at incidence. Second, results of experiments are presented in which the 
effects of boundary-layer interaction have been purposely suppressed by leaving 
a small gap between the body and the flare through which the boundary layer can 
escape. The resulting pressure distributions can then be compared directly 
with inviscid theories and thereby illustrate any deficiencies of the theories. 

Computer programs applicable to flared bodies at zero incidence have been 
developed at Ames Research Center (ref. 6) and elsewhere. A program based on 
the linearized characteristics method has also been used in references 7 and 8 
to calculate the first-order effects of angle of attack for smooth bodies. In 
the present report the linearized characteristics method is extended and 
applied to flows with embedded or secondary shocks. The results of sample cal- 
culations are compared with experimental pressure distributions obtained from 
the Ames 3- 5-Foot Hypersonic Wind Tunnel. Also, force coefficients and centers 
of pressure predicted by theory are compared with integrated experimental pres- 
sures and with other theories. 


SYMBOLS 


speed of sound 


A 

Ca 

% 

C P 

d 

h 

k b 

% 

K P 

Ke 

M 

P 

<1 




reference area, — : — 

k 

axial force coefficient (excluding viscous and base drag), 
axial force 

normal-force coefficient. 


normal force 


<3oo A 


p - p 

pressure coefficient, — ~ — 


diameter of cylindrical body segment 

enthalpy 

body curvature 

shock-wave curvature 

pressure gradient along a straight shock 
flow-angle gradient along a straight shock 
Mach number, ^ 
pressure 

pV 2 

dynamic pressure, 


2 


A -2222 


R radial distance to shock wave 

R*b blunt body nose radius 

S entropy 

s,n,t streamline coordinates (see fig, l) 

u velocity component in x direction 

v velocity component in r direction 

V scalar magnitude of velocity vector 

V velocity vector 

w velocity component in circumferential direction (crossflow velocity) 

W distance along a shock wave 

x™. distance to center of pressure 

x, r,$ cylindrical coordinates 

x,y, z rectangular coordinates 

X distance from center of spherical nose, X = x - 

a angle of attack 

(3 s/M 2 - 1 

7 specific heat ratio 

e 0 for two-dimensional flow, 1 for axisymmetric flow 

£ angle between shock wave and flow direction 

Q flow angle measured from x axis in meridional plane (see fig. l), 

tan ' 1 — 

| ,r| characteristic coordinates 

p density 

cr shock-wave angle measured from x axis 

cp crossflow angle (see fig. 1), sin -1 w/V 

$ azimuthal coordinate, cylindrical coordinates 
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Subscripts 


B body boundary condition 

d downstream shock conditions 

S shock boundary condition 

u upstream shock conditions 

o zero -order variable from solution of zero incidence flow field 

x first-order perturbation variable, implies a derivative with 

respect to a which is a function of x and r only 

00 free -stream conditions 

a derivative with respect to a 

Superscripts 

A,B, . . . points in the flow field 

1 coordinates fixed with respect to body axes 

tr coordinates fixed with respect to shock axes 


THEORY 


Linearized characteristics theory has been applied to a number of prob- 
lems involving smooth body shapes (see, e.g., refs. 7 to 9) . In this section 
the method will be applied to bodies with a compression corner which develops 
a shock embedded in a nonuniform flow field. The corner angle is restricted 
only by the condition that the shock must be attached and the flow behind the 
shock supersonic. The equations and boundary conditions are linearized with 
respect to the angle of attack, a. Thus the method yields information about 
the first-order changes of the flow field as a result of a small change in 
angle of attack. It depends on the prior solution of the flow at zero inci- 
dence and on prescribed initial conditions. In the present case the solutions 
for zero incidence are obtained with the computer programs described in refer- 
ences 6 and 10. Initial conditions for pointed bodies are obtained from cone 
theory as described in reference 8 . Initial conditions for spherically blunted 
bodies, to which present applications are restricted, depend only on the axi- 
symmetric solution which is calculated by the inverse method (ref. 10 ). 

The linearized characteristics equations will be outlined briefly in the 
following pages, and the boundary conditions for an embedded shock will be 
derived. The reader who does not care to follow the details of the derivation 
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may turn directly to the section on Calculation Procedure. One should, how- 
ever, take note of equations (l) below which define the notation used for the 
perturbation variables. A sketch of the coordinate axes and flow direction 
angles is presented in figure 1. 


Perturbation Equations 

The problem of calculating the perturbation flow field for a body with a 
compression corner differs from that for a smooth body only in the details of 
the embedded shock calculation. Therefore, before the embedded shock is dis- 
cussed the equations used in references 7 and 8 for smooth bodies will be 
listed. In the present form these equations apply to perfect or real gases in 
thermodynamic equilibrium. The subscript 0 is used to identify a quantity 
obtained from the zero angle -of- attack solution (a known quantity) , and the 
subscript 1 is used for the unknown perturbation variables defined by 


Pi(x,r) 


1 

COS $ 


p a ( x > r > 3 ’) 


1 

cos 0 



0i(x,r) 


cos $ 


0 a (x,r,«) 


l 

cos 0 





(i) 


cpi(x,r) = . - - cp„(x,r,$) = i 
^ * sxn 0 > ’ sin $ \oa/ a _ 0 

and so on. One should note that the perturbation variables are independent of 
the meridional angle, 0. Therefore the analysis can be made in any meridional 
plane; the leeward plane of symmetry is chosen in the present case. 


The basic perturbation equations, in terms of characteristic coordinates 
t> t\> are 


P 0 Sp 1 5© 

“iri - + = F - G 

p 0 v o 5r i 


(2a) 


where 


Pp 5p x _ 50! 

Po V o 2 31 " 


= E - G 


G = 


CPi + 0 COS 0 O 
Mor 


(2b) 
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Po 


MqPqV L 


, . Sp 0 SPo 

(A + 9±) ^ ~ (B + P o 0i) SnQ _ 


£1 | d 0 o _ d 
Mq \3so ° Sno 


E =. 


Po 
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3n, 


oj 
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and 


A = 


PoVo c 


Pi 


II + v ^ 3 (HI 


'|A| + WU 3 fK) + 5 o 2V ° 2 


IpA 


A/n 2 Pc 


B = 


Vo" 


5J |l + Cl 


'3h N 

.A/p 


V 

2P, 


Two additional equations give the gradients of the crossflow angle and entropy 
perturbations along the zero angle-of-attack streamline, s Q , which is the 
third characteristic direction. These are 


and 


dq>i _ ( 1 dPo _ sin g 0 \ Pi 

Sbq Vp 0 V 0 2 Ss 0 r J 1 rp 0 V 0 2 


dS Q 

as 0 _ " 01 dn 0 


(3) 


(4) 


The energy equation and the equation of state are needed to complete the system 
of equations. In terms of the perturbation variables these equations take the 
following form: 


Energy 


hi + VoVi = 0 


(5) 
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State 


Hi 




a i =■ 





Si - 



J 


( 6 ) 


The partial derivatives in equations (6) can be written explicitly for a per- 
fect gas but must be found numerically for a real gas as described in refer- 
ence 8. One supplementary equation that is needed for the present development 
is 

Vl = V Q Cpl (y) 

which gives the familiar crossflow velocity perturbation in terms of the cross- 
flow angle perturbation that appears in the previous equations. 

The initial conditions must specify the perturbation variables along a 
noncharacteristic line between the body and the shock in a supersonic flow 
region. At the present time, initial conditions are available only for pointed 
cones and spheres. In the present paper, application is made only to a spheri- 
cally tipped body, and the perturbation variables are all zero on the starting 
line in this case. 


Boundary Conditions for Smooth Bodies 

Although the equations given above are independent of the particular 
choice of reference axes, their boundary conditions are not. These conditions 
will vary depending on the choice of axes. In the present development the 
boundary conditions are given in terms of wind axes. Thus at the body one has 


( 0 i )]3 = -1 + (x 


R t>) 




(8) 


where 


(Si) B = K + [(x - Rb)eos 0 O + r sin 


0 ] 2§o 

0 ° J an 


K - 


0 

( cone 


for a spherical nose 
for a pointed nose 


(9) 
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At the bo w shock the perturbation conditions are 

(Pi)s = <*i - R i ^ 

/ \ R ^£o 

^s = 01 da Rl 8r 


( 01 )s = -x g - R i ^ 


(lOa) 

(10b) 

(10c) 



cos 0 O + cot Po sin 0c 


( 11 ) 


The two new perturbation parameters introduced in equations (10) and (ll) are 
the shock angular perturbation, a i, and radial perturbation, Ri. They are 
related geometrically by the relation 


cffii 

dx 


Ci sec 2 a 0 


( 12 ) 


The derivatives with respect to shock angle, a, can be obtained from the usual 
shock equations, and are evaluated in appendix A. 


Equations (l) through (12) completely specify the perturbation problem 
for smooth bodies. They apply also for cornered bodies except for the immedi- 
ate vicinity of the embedded shock. At a corner with an attached shock the 
term K in equation ( 9 ) will change, and all along the shock the jump condi- 
tions given in equations (10) and (ll) will be complicated by the nonuniform 
upstream conditions. These generalized shock conditions are developed next. 


p = p (r) 
p = p{r) 

9 = e[r) Shock wove for = 



Sketch (a) 


Boundary Conditions for an Embedded 
Shock 

Shock conditions f or Pjl,. Pi > <9_i. - 
Consider first a point on the shock 
some distance from the body and in the 
leeward plane of symmetry. The upstream 
conditions are nonuniform and the shock 
is displaced, as shown in sketch (a), 
as a result of pitching the body to 
angle of attack a. For the develop- 
ment of the perturbation shock condi- 
tions, it is convenient to introduce a 
derivative following the shock wave, 
that is, the rate of change of a quan- 
tity with respect to coordinates fixed 
with the shock wave. Thus the notation 
D/Da is used to indicate the deriva- 
tive 
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(13) 


Dp 

Da 


lim 

0(r>0 


P 8 


Po J 


a 


In terms of the original notation defined in equations (l), one may write, 
making use of a double prime superscript 

Pi"(x,r) 3 


<Pi"(x,r) = 


1 

Dp(x,r,$ ) 

cos $ 

Da 

1 

D6>(x,r,®) 

cos $ 

Da 

1 

Dcp(x,r,$) 

sin $ 

Da 


(14) 


The derivative D/Da can he written in a more suitable form by expanding flow" 
variables in a series in terms of the radial distance from point A. Using 
the pressure as an example, one has 


p 8 = p A + AR 


8p\ 

dr/ 


A 


+ 0(AR 2 ) 


(15) 


The perturbation expansions, in terms of wind axes, 

P A = P 0 A + ap x A + 0(a 2 ) 
AR =. aRi A + 0(a 2 ) 


( 16 ) 


are now introduced in equation ( 15 ), and the result substituted into equation 
( 13 ) • (Note that the subscripts 1 and a are synonymous here since the anal- 
ysis is made in the plane 0=0.) Terms of order a 2 are eliminated in the 
limiting process, and there results 


Dp -n 4. p ^ Po 

— = Pl + R X ( ^ 


(17) 


where it is understood that all quantities are evaluated at point A. This 
gives the first-order expansion of the shock derivative in terms of quantities 
evaluated in wind axes. Since the present calculations are made in terms of 
wind axes, the necessary shock conditions can now be easily formulated with 
the use of equation ( 17 ). The equations giving the jump conditions across a 
shock wave may be written generally as 


( 18 ) 


p d = 

p(Puv 

PHjv 

V u . 

£u) 

p d = 

p(Pu> 

Pu> 


£u) 

©u - 

e(Pu* 

Pu> 

V u , 

5u) 
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where 


- a - 0 U 


(19) 


and the subscripts d and u refer to conditions downstream and upstream from 
the shock -wave, respectively. Differentiating equations (l8) yields 


Dp d 

dp_ 

DPu 

+ 

^Pu 

+ 

DV u 

+ 


Da 

dp u 

Da 

SPu 

Da 

dv u 

Da 

^u 

Da 

p Pd _ 

dp_ 

DPu 

+ 

^Pu 

+ 

DV U 

+ Ifi- 

2k 

Da 

dPu 

Da 

dpu 

Da 

dv u 

Da 

d^ u 

Da 

D^u 

de 

D Ih 

+ 50 

^Pu 

be 

DV\, 

ds 

D^u 

Da 

^£*u 

Da 

dp u 

Da 

dv u 

Da 

d^u 

Da 


(20) 


These equations contain a k x 3 matrix of partial derivatives with respect to 
upstream conditions. Evaluation of these derivatives from equations (l8) 
is discussed in detail in appendix A. The presence in equations (20) of these 
upstream derivatives is one essential difference between the embedded shock 
conditions and equations (10) for the bow shock with uniform upstream condi- 
tions. In the case of uniform upstream conditions, only derivatives with 
respect to remain in equations (20), and these may be written, without 

loss of generality, as 


8^8 
d£u da 


( 21 ) 


since the partial derivative implies that 8 U is constant. However, all of 
the terms in equations (20) must he retained for the present problem. 

The last step needed to put the shock conditions into their desired form 
is to introduce equation (lj) into equations (20). The result (noting that 
boo/br - 0) is 


P 


id 



- Ri 



dp dPou _ dp dp QU ^ dp dV ou ^ bd ou ' \ 
dp u dr dp u dr dV u dr dp dr J 


dp 

P + 

1U dpu 


J iu 


d]2_ 

dpu 


+ V 


1U 


d£- - 

dV u 


dp 
1U da 


(22a) 
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djD _ _ / d p od ^ dp ou _ ^PpU dVpu Sp ^Qpi 

1<a - 1 ^ a 1 \ Sr Sp u dr dp u dr dV u dr da dr 


p ^- +D ^-+v §£- - o §£ 

P “ 3b,, Plu 3p u + T “ 3V a 0 . 


( 22 b) 


? id 


^ 5? - *i 


rp 

O 

CD 

ro 

he 

^Pou 

Se 

^Pou 

he 

^ v ou 

^Se h ^®ou 

L Sr 

5p u 

Sr 

hp u 

Sr 


3r + 

K^' 1 ) 3 rJ 


_ de_ 

Pl " 3P U 


4 . „ Ss_ . v de /’de , 
+ Plu d Pu + Vlu Sv u " e iuV§7 “ 1 


( 22 c) 


Equations (22) are the desired generalizations of equations (10). They contain 
additional terms which reflect the presence of upstream gradients in the basic 
zero incidence flow field, and upstream perturbations due to incidence. These 
upstream influences are multiplied by weighting factors determined from the 
shock equations. The weighting factors are the partial derivatives evaluated 
in appendix A. 

Shock condition s for the crossflow angle, cp^ .- The condition on the cross- 
flow angle is somewhat simpler due to the fact that cp = 0 everywhere for zero 
incidence. The present derivation is therefore basically the same as that 
given in greater detail in reference 7 . An additional term will appear, how- 
ever, which is due to the crossflow angle perturbation upstream of the shock. 

In terms of shock oriented coordinates, it is clear that 

Wxd = V lu ( 23 ) 

This equation states that the velocity component tangent to the shock is 
unchanged as required by the momentum equation. Use of equation ( 7 ) then 
gives 


Vou 

h'd = VT hu 

od 

The following equation, derived in reference 7 to transform cp 
wind axes. 


(24) 

from shock to 


9 i = 9 i" “ criXcos e Q - cot a Q sin e o ) 


(25) 
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may be applied to equation ( 2k ) with the result 


V, 


cpid + cri(cos 0 od - cot a 0 sin 0 od ) = [cp iu + 0x(cos 0 OU - cot cr 0 sin 0 OU ) ] 


V, 


od 


or 


Ud " 



cri 


cos 0 od - 


'ou 


V, 


od 


cos 0 OU ) - cot a 0 [ sin 0 od 


Vou 

Vod 


sin 0 



(26) 


For the bow shock wave 


and 


cpxu = 0 


'ou 


= 0 


so that equation (26) reduces to equation (ll). However, for a shock embedded 
in a nonuniform field, the more general expression (26) is needed. 


Shock conditions at t he corner . - The shock conditions derived above 
depend linearly on the shock angular perturbation, ox, and the radial pertur- 
bation, Rx. The calculation of these shock perturbation parameters for a 
general shock point will be outlined in the next section. Their initial 
values, at the surface of the body, are calculated in a different manner as 
outlined in the present section. 


Just as the shock conditions are most easily found with the use of shock 
oriented coordinates, the present analysis is simplified by body axes. Vari- 
ables expressed in body axes will be identified by a prime superscript. Refer- 
A. Position of corner for a = 0 ring to sketch (b) , one has 

B Position of corner for a>0 



Px' cos 


lim 

a-»o 



(27) 


In terms of body axes, the surface 
boundary condition is 

(6>i') b =0 (28) 

for both the upstream and downstream 
flow-angle perturbations. Making use 
of this condition, while differentiat- 
ing the shock relations (eqs. (l8)), 
one obtains 
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( 29 a) 


i dp t 3p dp * dp 

pia - ctx' ^ + v' ln + p* ^ + v lu 


da lu dpa P xu d^ 


dVi 


u 


t -4- ^e_ 4 . „» le- 


dp 


p id - 01 * 5a + p iu 5pu + p h 5^ + V “- Sv, 


U 


~ . 5© , 5© , 5© 

0ld = 0 = CTi s- + P, u 5— + p* ' 


+ v. 


50 


5a ^ m 1u ^ v iu sv, 


v u 


( 29 b) 

(29c) 


The partial derivatives used here are those introduced in equations (20) . 
Equation (29c) gives the initial value of the shock angular perturbation in 
terms of tody axes 


, -l {_% 5 0 . , 5e , „» 50 

ai = 50/5a V 1U ^ p iu 5p u + V i u 5Vi 


(30) 


The shock radial perturbation is, of course, zero when following the tody 

( R i f )b = 0 (3 1 ) 

Equations (29a) and (29h) together with the equation of state (eqs. (6)) yield 
the jump condition for the entropy perturbation 1 

^ (32) 


^d - I ~ 


dp/n Pld \dpA Pld 


P 

which, in turn, defines the new value for the constant, K, in equation (9) ; 
that is 


K = S 


id 


(33) 


downstream of the corner. 

Equations (29) to (33) complete the boundary conditions for the corner, 
but they must be converted to wind axes to be consistent with the rest of the 
analysis. This conversion has been discussed in detail in reference and 
only the results will be given here. These are: 


1 The condition, K = 0 in equation (9)j implies a constant entropy on the 
surface of a spherically tipped body at incidence. It is noteworthy that 
behind a corner the entropy varies as cos $ T in accordance with equation 
(32). This is an example of a situation where a body at incidence does not 
have a constant entropy everywhere on its surface. 
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c)p n 

Pi = Pi* + (x - Rb) 

Pi = Px* + < x - |p- 

. / \ 5 V 0 

Vi = V + (x - Rb) ^ 


- r 


- r 


dPp ^ 
Sx 

S Po 
Sx 

SVq 

a x 


R x = Ri* - R 0 tan a Q - (x - Rb) 

dob 


a x = p-l* - 1 - R c 


dx 


( 34 ) 


( 35 ) 

( 36 ) 


The transformations for entropy and flow angle have already been given in 
equations (8) and (9) above, and equation (26) for the crossflow angle applies 
here also without change. The derivative, dCo/^.? appearing in equation (36) 
is related to the shock curvature at zero incidence. It will not, in general, 
be zero when the upstream conditions are nonuniform. The derivative of the 
shock angle with respect to x can be obtained from the shock curvature which 
is evaluated in appendix A. 


This completes the specification of the perturbation shock conditions in 
a nonuniform stream. The next section presents a brief description of how the 
calculation of a typical shock point is performed. 


Calculation Procedure 

The equations and boundary conditions presented in the previous sections 
can be used to calculate numerically the entire perturbation field on a charac- 
teristic mesh evaluated for zero incidence. This calculation has been pro- 
grammed in conjunction with the axisymmetric program described in reference 6. 
After each point in the field is solved for a =■ 0 , the gradients of the axi- 
symmetric field are evaluated, and the perturbation equations solved. The cal- 
culation for smooth bodies is described in reference 7* In this section the 
details of the embedded shock calculation are outlined. 

Gradients of the a = 0 field. - The coefficients of the perturbation equa- 
tions and the boundary conditions contain derivatives, or gradients, of the 
zero incidence flow variables, p Q , p Q , 0 O , V Q , S Q . These gradients play an 
important part in the solution, and in the shock conditions in particular. For 
a general field point, a simple, linear backward difference method is suffi- 
ciently accurate to give good results. However, immediately behind a shock 
wave, such methods are not satisfactory, in general. A method which makes use 
of the shock curvature in conjunction with the equations of motion was there- 
fore adopted and is described in appendix A. 
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Corner point . - A typical charac- 
teristic mesh in the vicinity of a 
corner is shown in sketch (c). The 
fictitious point C (later to he dis- 
carded) is first calculated in the 
usual way by neglecting the corner. 
Conditions upstream of the shock at 
point B T are then calculated by 
linear interpolation between A-C, and 
at point F 1 by interpolation along 
D-E-C. Conditions at point B, down- 
stream of the shock, are then obtained 
by straightforward application of equa- 
tions ( 29 ) to (36). The calculation of 
conditions at point F is different from that for a general shock point since 
the field downstream of the shock is unknown at this stage in the calculation. 
An approximation to the shock angular perturbation, cti, at point F is 
obtained if the corner boundary condition (eq. (28)) is assumed to apply there. 
This assumes that dcxi/dx is small. Downstream conditions at point F can 
then be obtained with equations (29) to (36), and with corresponding upstream 
conditions. 



In the first attempts to perform the numerical calculations, the shock 
angle at point F was assumed equal to the value on the body at point B. 

This procedure, which neglects da 0 /dx, was used successfully in reference 6 
for zero incidence flows, but it caused large errors in the present perturba- 
tion solution. It was thus necessary to calculate the shock curvature at point 
B, and make use of this information to obtain a better approximation to the 
shock angle, a 0 , at point F. The equations developed in appendix A were used 
for this purpose. 


General shock point . - A typical characteristic mesh for a general shock 
point is shown in sketch (d). The fictitious point D is first calculated. 



and upstream conditions at point C* 
are obtained by linear interpolation 
along DE. Conditions at point A are 
also obtained by interpolation using 
points B and F (these are known from 
the previous shock point or, initially, 
from the corner point). One can now 
apply equation (2a) along the r\ coor- 
dinate from A to C, equation (l 2 ) 
along the shock from B to C, and equa- 
tions ( 22 a) and ( 22 c) at point C. 

These are written in finite difference 
form as follows: 

a(px C - Pi A ) + ( 0 i G - 0 i A ) = (F - G)At) 

( 37 ) 
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Eip = r x B + b(ai B + aiP) 


(38) 

(39) 


where 


Pi 


0i c 


_ C 

c 80 
0lC 


- Ri C d + e 

- Ri c f + g 


Po 

a = P0V0 2 

h = ^ (x^ - x®)sec 2 a Q 


m 


and d, e, f, and g can be identified by comparing equations (39) and (40) 
with equations (22a) and (22c). The four equations (37) to (40) in terms of 
the unknowns p^, OiP, cr-jC, and RjP can be solved for o^p* All of the 
remaining variables are then obtained by direct substitution. The procedure 
is iterated., making use of average values for the coefficients, until conver- 
gence of the pressure perturbation, Pi^, is obtained. 


The inviscid theory developed in the preceding pages is an exact numerical 
solution to the first-order angle -of- attack perturbation of the flow field; 
that is, it yields the exact initial slope of the pressure and other flow vari- 
ables. The largest angle of attack for which this linear approximation may be 
applied depends, of course, on the degree of nonlinearity of the actual flow. 
This nonlinearity can be assessed at present only by comparison with the proper 
experiment. An experimental program suited to this purpose has been conducted 
and will be described in the following section. 


EXPERIMENT 


In order to check the applicability of the perturbation theory, and to 
provide basic pressure data, a test program was conducted with the models shown 
in figure 2. A brief description of these tests follows. 


Facilty and Test Procedure 

The experimental study was conducted in the Ames 3. 5 -Foot Hypersonic Blow- 
down Wind Tunnel at Mach numbers of 5*2, 7*4, and 10. 5 and stagnation pres- 
sures of I3.6, 102.0, and 122.5 atmospheres (200, 1500, and 1800 psi), respec- 
tively, and at a constant stagnation temperature of 1166° K. The corresponding 
Reynolds numbers based on model length were approximately 0.9xl0 6 at M = 5.2, 
3.5xl0 6 at M =. 7-4, and 1.7xl0 6 at M = 10. 5. Data were taken at angles of 
attack of approximately 0°, 5°; and. 10 during an average total testing time 
of approximately 1 minute. A more general description of the test facility 
and its instrumentation is given in reference 11. A recent modification to 
the facility permits the test models to be inserted into the test section only 
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during the period when steady flow exists. This innovation was used for the 
present tests and greatly alleviates the effects of heating and stopping loads 
on the test models and instrumentation. 


Models and Apparatus 

The models were constructed of stainless steel and consisted of two inter- 
changeable noses, a cylindrical center section, and conical flare with a flare 
angle of 16.5 • The nose shapes were a hemisphere and a 22 - 1 / 2 ° half -angle 
cone. A photograph of the hemispherical-nose model is shown in figure 2(a) 
while a dimensional sketch, along with other details of the test models, is 
presented in figure 2(h)* The technique used in constructing the flare proved 
to he successful in bleeding off the boundary layer and, hence, in eliminating 
flare-induced separation on the cylindrical portion of the body (see ref. 12 ). 
It consisted of making the flare hollow and providing a gap between the flare 
and body cylinder. A gap size of 5 percent of the diameter of the body cylin- 
der was found to be suitable for the present tests. This gap size was slightly 
in excess of the maximum boundary-layer thickness at the flare estimated for 
the present models and test conditions. Preliminary tests were also made in 
which the gap size was varied from about 3 to 6 percent to determine if the 
pressures in the vicinity of the flare were affected. The results of these 
tests showed that the magnitude of the body pressures was not affected as long 
as the size of the flare gap was sufficiently large to prevent separation. 
Visible surface films were used in these tests to verify that boundary -layer 
separation had been eliminated with the chosen gap. 

The models were provided with thirty-nine 1.07 mm (0.042 in.) diameter 
pressure orifices distributed in two rows, 180 ° apart. Stainless -steel pres- 
sure tubing, 1.59 mm O.D. ( 0.0625 in.), was brazed into the walls of the 
models and led out through the base of the flare. About 3 feet aft of the 
flare, the stainless -steel tubing was joined to Viton high-temperature, flex- 
ible tubing, O .79 mm I.D. (O.O 375 in.) and 3.18 mm O.D. (0.125 in.), which led 
to the pressure cells. The total length of pressure tubing was about 4 meters. 
A table giving the location of the pressure orifices on the test models is 
presented in figure 2 (b). 


Instrumentation and Accuracy of Results 

The three ranges of strain-gage pressure cells employed in the tests were 
O. 34 , 0.68, and 3-^ atmospheres ( 5 , 10, and 50 psi). The cells were chosen to 
give maximum accuracy in measuring the pressures at a given model location for 
each test condition. Also, the pressure which was anticipated to be lowest 
was monitored on an oscillograph, and the data were recorded when this pressure 
stabilized. The precision of the final experimental data in coefficient form 
depends upon the accuracy in measuring the free-stream and orifice pressures 
and on the uncertainty in determining the free-stream dynamic pressure. The 
estimated over-all error in the pressure coefficient due to these uncertain- 
ties is ±0.007. The maximum estimated error in angle of attack, including 
possible -variations in free-stream angularity, was ± 0 . 3 °. 



RESULTS AND DISCUSSION 


Surface Pressures and Shock Shapes 

The experimental results for body surface pressures in coefficient form 
are presented in figures 3 and 4 for the hemispherical- and conical -nose 
models, respectively. Results for the hemispherical -nose model Q are o given for 
Mach numbers of 5.2, 7*^ and 10. 5 and at angles of attack of 0 , 5 , and 10 , 
while the results for the conical -nose model are given for Mach numbers of 5*2 
and 7.4 and angles of attack of 0 and 5 • 

Numerical calculations for the surface pressures on the spherical -nose 
model are also shown in figure 3. Theory and experiment are in agreement for 
a-0°, but differ, especially on the flare, for other angles of attack. The 
probable cause of this difference is discussed in the following paragraph. 

Since perturbation theory has not yet been developed for the expansion corner 
on the conical -nose model, the numerical results in figure 4 are only for 
a-0°, and agreement with experiment is excellent. In figure 5 calculated 
shock shapes are compared with experimental results from shadowgraph pictures 
for the blunt-nose model at a Mach number of 7*4. As the model is pitched to 
angles of attack of 5 and 10 > the bow shock remains unchanged relative to 
V^, while the flare shock must move with the body. The flare shock also moves 
relative to the body, becoming closer to the body on the windward side and 
farther removed on the leeward side. The calculated shock shapes on the lee- 
ward side and for a - 0 U agree well with experiment. On the windward side 
the experimental shock is somewhat steeper and farther from the body than the 
calculated result. 

Circumferential pressure dis tribution. - The comparison between theoretical 
and experimental flare pressures for the spherical -nose model at incidence 
(fig- 3) is not particularly good, even for 5 angle of attack. At 10 angle 
of attack the experimental pressure is considerably higher than predicted by 
the theory. On the other hand, the agreement on the spherical nose and on the 
forward section of the cylinder at a = 5 is very good. The reason for the 
disagreement between theory and experiment on the rearward portion of the body 
was studied in some detail and is attributed to the increased importance of 
second order (a 2 ) terms neglected in the theory. 

In order to illustrate this point, one must look at the second-order term 
in the expansion of the pressure. The series expansion in a and $ has the 
general form (see, e.g., ref. I3) 

p =. p Q 4- api cos 0 + a 2 (ps + P3 cos 2$) + 0(a 3 ) (4la) 


or 


p - p Q + api cos <& 4- a 2 [(p 2 - P3) + 2p 3 cos 2 0] + 0(a 3 ) (4lb) 

Thus if the pressure is plotted as a function of cos $ =- z, the second-order 
effect would appear as a displacement of the ordinate at 2 = 0 by a 2 (p 2 -p 3 ), 
and by a curvature (3 2 p/Sz 2 ) of 4a 2 p 3 . To determine if important second- order 
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effects were present in the data, additional tests were made with pressure 
orifices located in various 0 planes. These data are plotted versus cos 0 
in figure 6 for several x stations ahead of the flare; figure 6(a) is for 
a = 5° and figure 6(h) for a = 10°. Shown also are the linear pressure dis- 
tributions predicted by the theory. For 5 angle of attach the experimental 
variation is reasonably linear, although for x/d = 2.29 there is a tendency 
toward a nonlinear variation. For a = 10 , the first station (x/d - 0.45), 
which is on the sphere, shows a slightly nonlinear variation, and at rearward 
stations the experimental pressure distribution is increasingly nonlinear. 

There is a positive curvature to the pressure distribution at cos 0 = 0 which 
can be attributed to the term a 2 p 3 in equation (4Tb). In figure 7 similar 
data are shown for the flare pressures which also show this second-order effect, 
becoming very pronounced at a = 10°. On the flare these effects are amplified 
by the pressure rise across the shock. 

The discrepancy between theory and experiment is serious when detailed 
pressure distributions are desired. A need for more exact calculations, or 
estimates of the second-order term is clearly indicated. However, the linear 
theory may still be useful for predicting normal forces and moments at angles 
of attack much larger than that indicated by the comparisons with pressure 
data. The reason is that when equations (4l) are used for the integration of 
normal forces with respect to $, the second-order term is identically zero. 

Any nonlinearity in (% and (outside of viscous effects) must therefore be 
third order in a. These force coefficients and the center of pressure for the 
blunt -nose mode] are discussed next. 


Force Coefficients and Center of Pressure 

Integration of expe rimental pressures .- In order to evaluate the forces 
on the model at angle of attack it is, in general, necessary to integrate both 
in the circumferential and the longitudinal directions. For small angles of 
attack the circumferential pressure distribution has a cosine variation, and 
can therefore be integrated exactly. However, for the angles of attack of the 
present tests, it was found that the pressures deviate from the cosine curve. 
Nevertheless, by means of arguments given in the previous paragraph and illus- 
trated in figures 6 and 7; it is seen that the deviation is in the same direc- 
tion, and nearly equal, on the windward and leeward sides. In the case of the 
normal force and pitching moment, the additions to the cosine variation on 
windward and leeward sides of the body cancel. Therefore the experimental 
pressures in the plane of symmetry were used in conjunction with a cosine cir- 
cumferential distribution to calculate the forces and moments. A cosine dis- 
tribution was also assumed in calculating the axial force, even though the 
second-order effects on opposite sides of the body do not cancel in this case. 
As a result, the experimental data for should be slightly high at angle 

of attack. The error, however, is believed to be small since it occurs prima- 
rily on the flare which accounts for less than 10 percent of the total axial 
force at the test Mach numbers. Experimental data for the slope of the normal- 
force coefficient and the center of pressure at a = 0° were calculated with 
the 5° angle -of -attack data. 
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Comparison with the ory. - Aerodynamic coefficients obtained by integrating 
the experimental pressure distributions for the spherically blunted model are 
presented in figures 8 and 9* Figure 8 shows the effect of angle of attack 
for Mco =. 7.4, and figure 9 shows the variation with Mach number for a = 0 . 
Experimental results for the normal -force coefficient, axial-force coefficient, 
and center of pressure are compared with the present theory and with the theo- 
ries of references 2 and 3. The latter methods make use of a combination of 
simple theories. Jorgensen (ref. 2) assumes a linear pressure variation on 
the cylinder and applies cone theory to the flare. Seiff (ref. 3) makes use of 
blast-wave theory on the cylinder and Newtonian theory on the flare. Both of 
these methods as well as the present numerical results account for the nose 
contribution by means of modified Newtonian impact theory. 

In comparing results for varying angle of attack (fig. 8), one should note 
that the theories shown give only linear estimates. On the other hand, the 
experimental results exhibit the actual nonlinear nature of the aerodynamic 
parameters shown. This effect can be observed in the normal -force, axial -force, 
and center -of -pressure data which are all slightly high at 10 angle of attack 
relative to linear extrapolations of the lower angle data. Such trends are 
expected; however, attention is called to the previous section where reasons 
were given for expecting the experimental value for to be slightly high. 

For 5° angle of attack the experimental data are in best over -all agree- 
ment with the present numerical calculations. All the theories give essen- 
tially the same axial force, but there are sizable differences in normal force 
and center of pressure. The method of Jorgensen (ref. 2) predicts normal -force 
and center -of -pressure values about 10 percent higher than the present method; 
Seiff 1 s method (ref. 3) overpredicts the normal force by about 3O percent but 
gives the same center of pressure as the present method. However, the latter 
agreement is only for a Mach number of 7*4, as can be seen in figure 9* 

The comparisons at other Mach numbers between experiment and the various 
theories are shown in figure 9* Again good agreement is evidenced for the 
axial -force coefficient, owing mainly to the fact that the nose is the main 
contributor to this force. The present method gives the best agreement with 
experimental normal -force slopes, (^N a ) a _ 0 ^ and the estimates of reference 3 

are consistently high, while those of reference 2 tend to be slightly high at 
lower Mach numbers. With regard to the center of pressure, (x C p/d) a _ Q , th e 
situation is somewhat reversed. Seiff 1 s method gives the best over-all agree- 
ment with the experimental center of pressure, while Jorgensen 1 s method pre- 
dicts a curve nearly parallel to the data but about 8 percent too high. The 
present method predicts a greater forward shift of center of pressure with 
increasing Mach number than that indicated by experiment. 

As a consequence of overestimating the normal force, the method of ref- 
erence 3 also overpredicts the pitching moment. This result was observed in 
comparisons with data obtained from free -flight tests which were reported in 
references 4 and 14. The moment-curve slopes obtained from these free -flight 
tests decrease more rapidly, with increasing Mach number, than the predictions 
of Seiff* s theory (see ref. 14) • This trend is in general agreement with the 
center -of -pressure curve obtained from the present numerical results, but not 
with the trend of present experimental data, which shows nearly the same center 
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of pressure at Mach numbers of 7^ and 10.5* In view of this difference 
between present numerical calculations and experiment, additional data at 
higher Mach numbers are needed to verify the predicted shift of the center of 
pressure. 


CONCLUDING REMARKS 

The present study was undertaken to provide a better understanding of the 
inviscid flow over flared bodies at angle of attack. Pressure distributions on 
flared models, with hemispherical and conical noses, were measured at Mach num- 
bers of 5-2, 7.4, and 10-5* The flares on these models were provided with a 
means for boundary-layer removal in order to eliminate the large separated -flow 
region which would have been produced by the flare shock. The experimental 
pressure distributions therefore correspond to an idealized flow in which 
large-scale viscous effects have been suppressed. With the nature of the 
inviscid flow field thus established, the task of understanding flows compli- 
cated by viscous effects (which is the usual case) should be simplified. 

A linearized perturbation method was developed for the numerical calcula- 
tion of the flare shock conditions for small angles of attack, and the result- 
ing equations were incorporated into a computer program based on the linear- 
ized characteristics method. While this method can be regarded as being exact 
for sufficiently small angles of attack, comparison with the experimental 
pressure distributions showed significant deviations on the flare even for 5 
angle of attack. Since the experimental results were free of the effects of 
viscous separation, it was possible to show that the difference between theory 
and experiment was probably due to second-order (a 2 ) terms neglected in the 
present method. These second-order effects cause the circumferential pressure 
variation to depart from the cosine curve assumed in the perturbation theory. 
However, the departure from a cosine variation is an even function (i.e., 
cos 2 $) and therefore does not produce a cross force. (This is in contrast to 
the second -order effects studied in ref. 15 , which are caused by viscosity and 
which result in a cross force.) 

The mutual cancellation of pressures on windward and leeward surfaces 
thus results in normal forces and pitching moments which are more linear with 
a than would be anticipated from a knowledge of pressure distributions only 
in the plane of symmetry. The present theory was therefore useful for esti- 
mating the normal -force curve and the center of pressure for angles of attack 
of about 5°« Comparisons were made with more approximate theories, which are 
also linearized with respect to angle of attack. The present method was in 
best over-all agreement with the experimental results, but it predicted a 
greater decrease in stability with increasing Mach number than was measured 
experimentally. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Oct. 6, 1965 
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APPENDIX A 


EVALUATION OF SHOCK DERIVATIVES AND FLOW GRADIENTS 
BEHIND A SHOCK 


Calculation of the perturbation flow over a body at incidence requires the 
knowledge of gradients of the zero incidence field* Evaluation of the field 
gradients immediately behind a shock wave is accomplished with the use of the 
equations of motion in conjunction with certain partial derivatives which were 
introduced in the shock boundary conditions. These partial derivatives 
describe the downstream influence of unit changes in the shock angle and condi- 
tions upstream of the shock. The calculation of the shock derivatives from 
the oblique shock equations is first described. The equations for flow-field 
gradients are then derived and their application to the calculation of the 
shock -wave curvature is illustrated. 

SHOCK-WAVE DERIVATIVES 


It is convenient to use index notation where f^ are the unknown condi- 
tions downstream of the shock wave, and g^ are the known upstream conditions. 
Then the shock equations which determine the downstream conditions may be writ- 
ten symbolically as (cf. eqs . (l8)) 

fi = fl(gl,g2.>g3,g4) (Al) 

In particular, the variables are identified as 

fi = p,p ,9 for i = l-*3 

811(1 g k = Pu.Pu^u^u f °r k = 


where 


m ° 


Differentiating (Al) with respect to a, as in equations (20), gives 

D£i ^ y Sfi Dgk 
Da ~ dg k Da 
k=i 


(A2) 


where the tensor df i/Bgfc consists of the 4x3 matrix 
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(A 3 ) 


equations given in reference 1 6 are used to obtain explicit relations for the 
elements of (A3). These are listed below: 


5 p _ f 7 - 1\ 

^Pa = " V + V 


3 p 

5 Pu 


27/(7 + 1) 

V u 2 sin 2 a 

p u 7 

\7 + V V + V p u V u 2 Sin 2 0 . 


3 e 

bVu 


-sin 2 


G tan a 


(1 + l\ Pu v u 2 

v 27 y p u 2 


D V «= 
M uu 

7Pu 


sin 



3 p 2 V u g sin 2 a 
b^ = 7 + 1 


(A4) 
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(A8) 
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p u 

1 

K, 

PuV u 2 J 


(Al4) 


be 


dg sin(cr - 0 )cos(cr - d) 4/Pu cos 2 (a - 0 ) 

da - sin a cos a t . \ ,t 2 -, 'S- 


(7 + l)p u V u sin a 


(A 15 ) 


For a real gas in equilibrium, the shock conditions are implicit (i.e., 
they must be obtained by iteration), therefore the derivatives above must be 
evaluated numerically. Although the results of flow-field calculations pre- 
sented herein are for thermally perfect air, it is instructive to compare equa- 
tions (a 4) to (A12) with the results of numerical calculations for dissociated 
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air- This comparison is made in figure 10 as a function of the upstream veloc- 
ity normal to the shock. The gas properties used in the calculations were 
obtained from curve fits to equilibrium thermodynamic properties generated in 
the manner described in reference 17- 


Gradients of the Zero Incidence Flow Field 

The derivatives of flow variables along streamlines and normal to the 
streamlines play an important part in the calculation of the perturbation flow 
field. For the usual mesh point, a simple linear backward difference method 
is used for evaluating these derivatives. This method is not well suited, 
however, for calculating the downstream gradients along the bow shock, and 
especially along an embedded shock. For these cases, a more accurate method 
is used which makes use of the differential equations of motion directly. 

These are written in terms of streamline coordinates (see, e.g., Hayes and 
Prob stein, ref. 18) 


g 2 dp bo sin Q 

pV 2 5s ~ " Bn 6 r 


(Al6) 


where e = 0 for 
e = 1 for 


l 5p = 

pV 2 5n 5s 

two-dimensional flow 

axisymmetric three-dimensional flow 


(A17) 


and s, n are distances along and normal to the streamline. The vector 
gradients of p and 0 along the shock wave may be written in terms of com- 
ponents along the s, n coordinates in the following way (see sketch (e)) 


n 



Sketch (e) 


dp 

dW 


dp 

cos £ — + sin t 
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dp 

dn 


(A18) 


d0 

dW 


d 0 , . ,, 

cos £ + sm £ 


de 

dn 


(A19) 


Equations (Al6) to (A19) contain four 
unknown derivatives of p and Q with 
respect to s and n in terms of 
gradients along the shock. Substitu- 
tion yields the following expressions 
for the streamline pressure gradient 
and the streamline curvature : 


la. _ 
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cos £ + pV 2 sin t, 



sin 0 sin £ 
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cos 2 £ - [3 2 sin 2 £ 


(A20) 


5© 


b 2 ^ . 

— ^ r — sin C + COS L 

pV 2 5w b _ b 



sin 0 sin £ 


5s 


cos 2 t, - |3 2 sin 2 £ 


(A21) 


The remaining derivatives can then be found from equations (Al 6 ) and (A17). 
The entropy gradient can also be expressed in terms of the gradient along the 
shock. In this case, since the entropy is constant on streamlines, one 
obtains 


|S = _1 H (A22) 

dn sin £ &W 

Equations (A20) and (A2l) require derivatives of p, 0, and p along the 
shock wave, evaluated on the downstream side. These can be obtained in terms 
of the upstream gradients, which are presumed known, by means of the shock 
derivatives (A 3 ); that is 


5p 

y S P Sg k 

5w = 

L 5g k Sw 


k=i 


4 

5e 

y 50 s gk 

5w 

Zj 5g k 5w 


k=i 


(A23) 


(A24) 


where the derivatives 5g k /5W are determined from the upstream gradients. 

Equation (A23) contains the shock curvature, 5a/5w, which enters through 
the derivative of £ u 


^£u 5a 
5w - 5w 5w 


For a general shock point the shock curvature is found by a linear difference. 
Near a corner, however, this quantity is evaluated differently as described 
next. 
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Shock-Wave Curvature 


When a nonuniform supersonic stream flows over a compression corner, the 
initial shock angle at the corner is equal to the two-dimensional oblique 
shock value, and the shock will have a finite curvature even though the body 
may be straight. On the other hand, a curved body will produce a curved shock 
in a uniform stream. The latter problem has been considered by a number of 
authors (see, e.g., ref. 19)- In either case, the initial shock curvature can 
be evaluated from the equations given in this appendix. The shock curvature 
depends on the streamline curvature downstream of the shock which is determined 
from a knowledge of the body shape. Equations (A 23 ) and (A24) can be used in 
conjunction with equation (A2l) for the streamline curvature in evaluating the 
shock-wave curvature. First, however, it will be convenient to rewrite equa- 
tions (A 23 ) and (A24) in the form 


dp dp 

§w = + ^ % 


(A25) 
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dw 
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and where 


K W = 


bo 

b W 


(A26) 


is the shock -wave curvature. The terms Kp and K q represent the gradients 
along a straight shock wave, and the remaining term represents the gradient 
due to shock curvature. Substitution of equations (A25) and (A2 6 ) into (A2l) 
yields, after some rearranging 


Kp(cos 2 £ - (3 2 sin 2 £) - cos £ 


% = 




\pWj K P tan ^ + K 0 + r sin Q sin 


p 2 . ,, 3s ,, 

z—s — sin t + cos c, 
p V 2 bo 00 


(A2J) 


where 
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is the body surface curvature "behind the shock. This expression depends on 
the body curvature, on the flow-field gradients upstream of the shock, and on 
the derivatives of the shock-wave equations. These conditions are known from 
the solution of the flow upstream of the corner. The initial shock curvature 
for a two-dimensional or axisymmetric flare in a nonuniform stream is there- 
fore determined. 
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(a) Photograph of hemispherical -nose model mounted in 3« 5-foot HWT test section. 

Figure 2.- Wind-tunnel models. 
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(b) Dimensions. 
Figure 2.- Concluded. 



(a) Moo = 5-2 

Figure 3*“ Surface pressures for the hemisphere -cylinder -flare model 
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Figure 6.- Concluded. 
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Figure Comparison of experimental and theoretical estimates of force 

coefficients for a hemi sphere -cylinder -flare ; a = 0 . 
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(a) Derivatives with respect to upstream pressure 
Figure 10.- Shock wave derivatives. 
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